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摘 要 :为 进一步 验证 所 提 人 饱和 砂 土 三 前 弹 塑 性 边界 面 本 构 模 型 的 适用 性 ,以 南昌 地 区 饱和 砂 土 为 
研究 对 象 ,利用 PFC3D 颗粒 流 建立 了 单调 及 循环 荷载 作用 下 不 同 含 泥 量 饱和 砂 士 外 训 柔 性 膜 的 党 
规 三 轴 和 真 三 轴 数 值 试 验 模 型 。 常 规 三 轴 条 件 下 的 室内 试验 .PFC3D 模拟 及 理论 模型 间 的 结果 对 
比 验 证 了 3 者 之 间 在 固 结 排 水 时 土 的 偏 应 力 与 轴 向 应 变 、 轴 向 应 变 与 体 应 变 , 以 及 固 结 不 排水 时 土 
的 偏 应 力 与 轴 向 应 变 、 轴 向 应 变 与 孔隙 水 压力 之 间 的 关系 均 具有 较 好 的 一 致 性 ,达到 了 可 以 进行 数 
值 试验 的 程度 ,表明 所 建 PFC3D 模型 可 以 代替 各 种 复杂 加 载 条 件 下 的 饱和 砂 土 室内 试验 。 另 外 ， 
将 PFC3D 真 三 轴 固 结 不 排水 数值 试验 和 理论 模型 间 的 结果 作 了 对 比 ,研究 了 不 同 含 沁 量 饱和 砂 土 
的 中 间 主 应 力 影响 问题 。 其 中 ,在 静 真 三 轴 固 结 不 排水 条 件 下 ,单调 荷载 下 的 PFC3D 模拟 与 中 间 
主 应 力 影 响 系 数 =0.5 时 的 理论 曲线 整体 吻合 均 较 好 ,其 曲线 误差 在 轴 向 变形 前 期 基本 保持 5% 
以 内 ;在 后 续 的 动 真 三 轴 边 界面 理论 模型 验证 时 ,循环 荷载 下 的 理论 模型 计算 均 采 用 了 2 =0.5。 在 
动 真 三 轴 固 结 不 排水 条 件 下 ,PFC3D 数值 模拟 和 理论 计算 结果 在 应 力 - 应 变 、 循 环 次 数 -孔隙 水 压力 
关系 曲线 也 比较 吻合 , 含 泥 量 越 多 应 变 发 展 速度 越 快 、 动 剪 切 模 量 更 小 、 沿 回转 面积 更 大 , 孔 队 水 压 
的 振动 幅 值 和 绝对 值 都 随 之 变 大 。 
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Approaches to the triple-shear elasto-plasticity bounding 
surface model for saturated sands (3) 
—PFC3D numerical test verifications 


HU Xiaorong, CAI Xiaofeng,LIU Cao 


(School of Infrastructure Engineering, Nanchang University ,330031 Nanchang China) 


Abstract:In order to verify the applicability of the triple-shear elastoplastic bounding surface constitutive 
model for saturated sands proposed previously, the conventional triaxial and true triaxial numerical test 
models using PFC3D were established for Nanchang saturated sands with different clayey contents wrapped 


in flexible membranes and bearing the monotonic and cyclic loading. Comparisons among the results from 
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the laboratory conventional triaxial CD and CU tests, the PFC3D simulations, and the theoretical models 
showed good consistency in relationships existing in the deviatoric stress vs the axial strain and the axial 
strain vs the bulk strain under the CD conditions , and the same results existed in the deviatoric stress vs 
the axial strain and the axial strain vs the pore water pressure under the CU conditions. The conclusions 
mentioned above indicate that the numerical simulations with built PFC3D models can be used to verify the 
theoretical models replacing the laboratory tests for saturated sands under various complex loading condi- 
tions. In spite of that, results between the PFC3D true triaxial numerical tests and the theoretical models 
are compared to study the effects of the intermediate principal stresses for the saturated sands with different 
clayey contents. Because the curves deduced from the PFC3D simulation acting on the monotonic loads co- 
incide with those from the theoretical models in the static true triaxial CU tests when b 20. 5 , giving the er- 
ror within 596 ,b 20.5 is adopted for all to cyclic load computations for theoretical models in the later true 
dynamic triaxial verifications to the theoretical models of bounding surface. Result comparisons between the 
PFC3D simulations and the theoretical models show good curve consistency in the stress vs the strain and 
the cycle times vs the pore water pressure under the dynamic true triaxial CU conditions. The greater clay- 
ey contents lead to the faster strain developing speed ; the smaller dynamic shear moduli , the larger hystere- 
sis loop areas,the greater vibration amplitudes and the absolute pore water pressures. 

Key words:saturated sand ; triple-shear elasto-plasticity bounding surface constitutive model; PFC3D nu- 


merical test verification ; monotonic load;cyclic load 


南昌 地 铁 开 挫 并 运行 在 他 和 砂 土 层 中 ,为 了 研 
究 其 在 动静 价 载 作用 下 的 力学 特性 ,文献 [12] 提 出 
局 忠 和 砂 土 的 三 芒 弹 塑性 边界 面 本 构 模 型 并 做 了 背 
规 动静 三 轴 试 验 验 证 ,结果 表明 理论 模型 与 常规 三 
轴 汇 验 结 末 较 吻合 ;遗憾 的 是 ,未 能 对 理论 模型 作 相 
RHE 、 动 真 三 轴 试 验 验 证 。 然 而 ,该 试验 可 能 会 因 
议 但 价格 昂贵 而 难以 实现 ; 妃 外 ,试验 工作 量 大 、 技 
术 蔷 求 高。 数值 试验 能 从 多 角度 去 揭示 材料 的 力学 
特 和 性 及 破坏 机 理 , 可 以 代 蔡 各 种 复杂 甚至 是 无 法 实现 
加 载 条 件 的 室内 试验 。 如 有 果 能 合理 地 选用 软件 .正确 
地 建立 计算 模型 ,准确 地 确定 模型 计算 参数 并 经 得 起 
与 典型 真实 试验 及 理论 模型 的 验证 , 则 可 认为 该 计算 
软件 下 的 模型 达到 了 可 进行 数值 试验 程度 ,并 接近 真 
实 试验 ,可 以 代替 真实 试验 来 验证 理论 模型 。 

土 体 是 由 离散 土 颗粒 为 骨架 组 成 的 多 相 体 , 基 
于 土 颗 粒 尺 度 的 三 维 颗粒 流 软件 PFC3D 在 材料 力 
学 特性 数值 模拟 方面 得 到 广泛 应 用 。 文 献 [9] 
分 析 了 人 砂 土 的 破 肆 特性 对 宏观 力学 特性 的 影响 ; 文 
献 [10-13] 利 用 PFC2D 分 析 了 笑 土 和 砂 土 的 细 观 力 
学 特性 及 破坏 过 程 ;文献 [14] 人 研究 了 微 重力 作用 下 
振动 容 副 中 颗粒 的 输送 问题 ;文献 15] 探讨 了 无 盘 
桩 路 堤 中 拱 的 谷 载 恢复 机 制 , 进 行 了 数值 和 物理 样 
本 之 间 的 效率 啊 应 ,体积 变形 和 位 移 特 性 的 比较 ; 文 
献 L16 研究 了 溶 胀 颗粒 状 材料 的 孔 隐 度 、 应力 状 态 


和 溶 胀 度 之 间 的 关系 ;文献 [17] 曾 采用 PFC3D 对 南 
昌 地 区 饱和 砂 土 做 了 单调 伍 载 下 的 三 轴 试 验 模 拟 并 
与 实际 试验 结果 作 了 对 比 ,但 建 模 时 采用 了 刚性 墙 
去 模拟 真实 试验 中 的 橡皮 膜 , 这 与 实际 情况 不 符 。 

为 进一步 验证 文献 [|1] 所 提 人 饱和 砂 土 三 羡 弹 逆 
性 边界 面 模 型 的 适用 性 ,在 文献 [1-2,17 |] 基础 上 对 
南昌 地 区 饱和 砂 土 进行 PFC3D 数值 试验 验证 。 数 
值 试 验 为 单调 衔 载 作用 下 不 同 含 泥 量 饱和 砂 土 外 理 
柔性 膜 的 常规 三 轴 固 结 排水 和 固 结 不 排水 试验 和 真 
三 轴 固 结 排 水 试验 ,循环 谷 载 作 用 下 动 三 轴 和 动 真 
三 轴 的 固 结 不 排水 试验 。 研 究 特 点 是 :中 数值 试验 
中 的 饱和 人 砂 土 外 于 柔性 膜 模 型 能 够 模拟 真实 试验 中 
的 橡皮 膜 ;@ 思 对 文献 [1 中 的 理论 模型 可 进行 较 完 
整 的 动 . 毅 真 三 轴 数 值 试 验 验证 。 具 体内 容 如 下 。 

1)PFC3D 数值 模型 建立 及 细 观 参数 标定 。 采 
用 并 编写 了 PFC3D 柔性 外 膜 程序 来 更 好 地 模拟 真 
实 三 轴 试 验 中 土 样 的 外 右 橡 皮膜 ,可 以 有 效 避 免 原 
有 PFC3D 数值 模型 中 刚性 墙 所 导致 土 样 上 下 变形 
相同 的 缺点 。 数 值 模型 中 土 样 的 细 观 参数 以 文献 
[2 的 一 部 分 室内 三 轴 试 验 结果 确定 。 

2) 饱和 人 砂 土 常规 三 轴 试 验 结 采 对 比 。 利 用 
PFC3D 数值 模型 来 模拟 文献 [2] 的 其 他 室内 常规 三 
轴 试 验 , 并 同时 与 文献 [1j] 中 的 理论 计算 结果 进行 
三 者 对 比 验 证 ,以 验证 PFC3D 建 模 及 其 土 性 细 观 参 
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数 的 可 行 性 。 

3) 饱 和 砂 土 真 三 轴 结 采 对 比 。 采 用 内 容 2) 验 
证 了 PFC3D 数值 模型 和 土 性 细 观 参数 进行 动静 真 
三 轴 模 拟 , 并 与 理论 模型 计算 结果 进行 对 比 ,以 进 一 
步 验证 文献 [1 所 提 本 构 模型 的 正确 性 。 


1 PFC3D 模型 建立 及 细 观 参数 标定 


1) 基 于 文献 [2] 室 内 三 轴 试 验 的 实际 砂 土 情况 
和 颗粒 得 分 试验 结果 。 首 先生 成 超过 5 000 个 数值 
砂 土 颗粒 , 且 砂 土 颗粒 半径 为 实际 砂 土 样 的 3. 5 
倍 " ”1。 数 值 砂 土 的 初始 孔隙 率 与 室内 三 轴 试 验 一 
致 , 均 为 0.3 ,依照 饱和 砂 土 含 泥 量 的 不 同 (0% 、 
546.1090 ) 来 设置 砂 土 样 体积 占 比 ,试验 所 采用 的 饮 
和 和 束 土 与 数值 砂 土 的 级 配 曲线 见 图 1。 其 中 ,X 为 小 

粒 径 土质 量 百 分 数 ,d 为 粒 径 。 根 据 不 同 排水 
条 件 ,分 别 进行 固 结 排水 试验 (CD 试验 ) 和 固 结 不 
íi CU 试验 ) 7" 。 针 对 试验 过 程 中 可 能 出 现 
的 种 型 加 载 速度 突变 而 导致 的 因数 值 砂 土 受 力 不 均 

坏 的 现象 ,采用 的 方法 是 :在 柔性 外 膜 内 通过 缩 
旋 沙 土 颗粒 的 大 小 以 控制 砂 土 颗粒 不 外 溢 , 并 且 在 
FOREZ SEPA 10 次 为 一 循环 对 砂 士 颗粒 速 
度 进行 清 零 , 从 而 确保 计算 稳定 。 
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图 1 饱和 砂 土 与 数值 砂 土 的 级 配 曲 线 
Fig.1 Grading curves of saturated sands and 
Ee A ee ge 
2 ) 饱 和 人 砂 土 的 常规 三 轴 试 验 数值 模型 建立 。 如 
图 2 所 示 , 设 置 一 个 78.8 mm x 78. 8 mm x 160 mm 
的 长 方 体 数值 模型 边界 空间 ,根据 文献 [2 ] 的 饱和 
人 砂 土 试 样 实际 尺 寸 在 边界 空间 内 建立 直径 为 
39. 1 mm, 高 为 80 mm 的 圆柱 体 土 样 。 
3) 砂 士 真 三 轴 数 值 模型 的 建立 。 将 常规 三 轴 圆 
柱 模型 改 为 六 面体 模型 ,模型 尺寸 为 80 mm x 
39.1 mm x39. 1 mm, 如 图 3 所 示 , 其 他 设置 与 常规 
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(a) 模型 边界 


(b) 含 泥 量 0% 


(c) EES% 
图 2 常规 三 轴 数 值 砂 土 试 样 模型 


Fig.2 Conventional triaxial sample models 


(a) 真 三 轴 边 界 


(d) 含 泥 量 10% 


(b) 含 泥 量 0% 


(c) 含 泥 量 5% (d) 含 泥 量 10% 
图 3 真 三 轴 试 样 模型 
Fig.3 True triaxial sample models 
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PFC3D 数值 模型 会 导致 砂 土 在 受 力 后 发 生 沿 侧 回 的 
均等 变形 ,不 符合 常规 三 轴 试 验 中 土 样 爱 外 右 橡 皮 
膜 约束 时 的 柔性 变形 特征 。 为 了 保证 PFC3D 数值 
村 型 能 够 反映 试验 中 橡皮 膜 的 柔性 加 压 功 能 ,开发 
了 PFC3D 柔性 外 膜 程序 。 其 中 , 固 结 排水 试验 中 的 


柔性 外 膜 模拟 是 以 颗粒 粒 径 为 现实 试验 人 砂 土 最 小 颗 
粒 粒 径 的 2/3 来 生成 条 性 颗粒 膜 ,采用 接触 暮 结 模 
型 ,以 柔性 外 膜 颗 粒 的 初始 位 置 不 断 循环 来 创建 并 
排列 好 彼此 相 切 的 颗粒 ,如 图 4(a) ~ 图 4(c) 所 示 。 
回 结 不 排水 试验 中 的 柔性 外 膜 则 是 在 生成 颗粒 后 在 
刚性 墙 的 基础 上 内 车 一 层 柔 性 外 膜 ,两 者 共同 组 成 
数值 砂 土 样 的 边界 面 ,为 了 克服 刚性 墙 无 法 移入 接 
触 这 一 问题 ,在 加 压 过 程 中 并 不 是 以 柔性 外 膜 桥 粒 
给 予 等 效 集中 力 ,而 是 通过 刚性 墙 来 控制 围 压 , 如 图 
4(\d) 和 图 4(e) 所 未。 
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N 图 4 柔性 外 膜 的 建立 过 程 


> Fig. 4 Establishment process of flexible outer film 

S25) 饱 和 砂 土 细 观 参数 的 确定 。 以 试验 围 压 
200kPa A JEE 5% 的 饱和 砂 土 的 循环 动 三 轴 试 验 
( 取 振 动 次 数 的 前 20 次 ) 为 对 象 进行 细 观 参数 标定 ， 
通 竹 对 砂 土 的 接触 刚度 敏感 性 分 析 (2 x 107 N/m, 
1 X10" N/m 和 5 x 105 N 和 m) 以 及 柔性 外 膜 颗粒 的 接 
触 刚度 敏感 性 分 析 (5 x10 N/m,1 x10* N/m 和 2 x 
10* NMm) 可 知 ,采用 泥 颗 粒 刚度 1 x10% N/m 、 砂 颗粒 
刚度 1 x10 N/m 以 及 柔性 外 膜 颗 粒 刚度 2 x 10 N/m 
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时 与 试验 结果 最 为 匹配 。 由 数值 模型 的 阻尼 人 敏感 性 
分 析 (0.4.0.7 和 0.9) 可 知 , 阻 尼 设 置 为 0.7 时 更 为 
恰当 。 经 过 不 断 调整 计算 参数 ,饱和 砂 土 样 的 最 终 
PFC3D 数值 模拟 细 观 参数 见 表 1 。 

利用 表 1 所 确定 的 饱和 砂 士 细 观 参数 来 模拟 试 
验 围 压 100 kPa 和 含 泥 量 5% 的 饱和 人 砂 土 固 结 排水 
数值 试验 以 及 试验 围 压 200 kPa 和 含 泥 量 5% 的 饱 
和 砂 土 固 结 不 排水 循环 动 三 轴 试 验 ,PFC3D 的 模拟 
结果 与 静 动 三 轴 试 验 对 比 结果 如 图 5 所 示 。 
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(a) CD 下 偏 应 力 g- 轴 向 应 变 g, (b) CD 下 体 应 变 s.- 轴 问 应 变 8， 


一 试验 结果 
go- 一 PFC 模拟 结果 


一 试验 结果 
一 PFC 模拟 结果 


(c) 循环 荷载 下 的 偏 应 力 g- 轴 癌 应 变 sg 


(d) 循环 荷载 下 的 孔 院 
IK Eu- AKN 


(e) CD 试验 数值 模拟 压缩 结果 (室内 CD 试验 压缩 结果 
图 5 试验 结果 与 模拟 结果 对 比 


Fig.5 Comparisons between tests data and simulation results 


表 1 砂 土 模型 细 观 接触 参数 汇总 


Tab.1 Summary of meso contact parameters of sand models 
"— "— odd pad icd ipod mm 阳 尼 us m 
砂 线性 接触 1 x10” 1 x10” - - 0.70 0.7 2 600 
泥 粒 Bof iih ds 1 x 10^ 1 x10° 331097 3x10 0.24 E 2710 
IRINA Bof ia 1 x 10^ 1 x 10^ i999 iod 0.70 0.7 1 000 
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由 图 5(a) 和 图 5(b) 可知, 固 结 排水 三 轴 试 验 
结果 与 PFC3D 数值 模型 模拟 结 末 均 较为 接近 ,数值 
模型 能 较 好 地 反映 饱和 含 泥 砂 土 的 软化 阶段 以 及 体 


积 收缩 现象 。 由 图 $(c) 和 图 $(d) 可 知 ,PFC3D 数 
值 模型 能 够 很 好 地 模拟 固 结 不 排水 循环 动 三 轴 试 验 
Hans [e] PS ,累积 塑性 应 变 和 和 孔 辽 水 压力 随 循环 次 数 
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的 变化 也 较为 一 致 。 由 图 5(e) 和 图 51 可 知 , 在 回 
结 排水 三 轴 试 验 中 ,PFC3D 数值 模拟 所 得 土 样 剪 切 
破坏 形状 与 实际 试验 试 样 是 相似 的 。 上 述 结 末 表 明 
本 研究 所 建 模型 及 所 标定 的 模型 计算 参数 是 正确 
的 ,可 用 于 对 饱和 人 砂 土 作 进一步 的 力 竺 特性 分 析 。 


2 饱和 砂 土 常规 三 轴 CD 和 CU 试验 


结果 对 比 


为 了 进一步 验证 本 人 研究 所 建 PFC3D 模型 的 正 
确 性 ,并 将 其 作为 数值 试验 机 来 验证 文献 [1] 中 理 
论 模 型 的 可 能 性 ,用 所 建 PFC3D 数值 模型 对 围 压 分 
别 为 100 kPa、200 kPa 和 400 kPa、 含 泥 量 分 别 为 


0% 5% 和 10% 的 饱和 人 砂 土 做 了 固 结 排 水 (CD) 与 
结 不 排水 (CU ) 试验 模拟 ,并 与 室内 试验 和 理论 计算 
结果 做 了 对 比 。 模 型 采用 圆柱 体 , 通 过 用 Fish 语言 
代码 编程 实现 如 下 加 载 :首先 在 相应 的 排水 条 件 下 ， 
对 PFC3D 数值 土 样 的 3 个 方向 同时 加 载 ,直至 目标 
最 小 主 应 力 ; 符 完 成 固 结 后 ,保持 目标 最 小 主 应 力 大 
小 不 变 , 在 紧 直 方向 分 次 加 载 至 目标 最 大 主 应 力 , 直 
至 数值 试 样 剪 切 破 坏 。 其 中 , 固 绪 不 排水 模拟 试验 
中 的 侧身 围 压 增 量 视 为 实际 试验 中 的 筷 际 水 压 。 图 
6 和 图 7 为 固 结 排水 试验 中 轴 回 应 变 分 别 与 俩 应 
力 \ 体 应 变 的 关系 曲线 ,图 8 和 图 9 为 固 结 不 排水 试 
验 中 轴 向 应 变 分 别 与 仿 应 力 、 筷 际 水 压力 的 天 系 
曲线 。 


-a- p-100 kPa 下 试验 曲线 -»- p-100 kPa 下 试验 曲线 
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试验 与 理论 和 实际 试验 结果 非 第 接近 ,能 够 较 好 地 
反映 不 同 含 泥 量 饱 和 砂 士 的 软化 特性 ,并 与 文献 
[21 ] 的 绪 采 相同 。 由 图 7 可 以 看 出 , 固 结 排水 模拟 


图 7 CD 试验 的 体 应 变 与 轴 疝 应 变 关 系 曲线 
Fig.7 Curves of the bulk strain vs the axial strain in CD tests 
由 图 6 可 以 看 以 看 出 ,PFC3D 的 固 结 排水 模拟 。 试验 中 的 体 应 变 模 拟 曲 线 与 实际 试验 结果 和 理论 计 
算 结 朱 也 基本 符合 ;对 于 不 同 含 泥 量 和 围 压 下 的 饱 
和 砂 土 ,其 受 力 后 的 体积 收缩 和 体积 膨胀 现象 也 能 
得 到 较 好 地 体现 。 


-a- p-100 kPa 下 试验 曲线 
-e- p-100 kPa 下 PFC 模拟 曲线 


-a- p-100 kPa 下 试验 曲线 
-e- p-100 kPa 下 PFC 模拟 曲线 
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æ- p-100 kPa 下 试验 曲线 
-e- p-100 kPa 下 PFC 模拟 曲线 
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图 8 CU 试验 的 侦 应 力 与 轴 癌 应 变 关系 曲线 


Fig.8 Curves of the deviatoric stress vs the axial strain in CU tests 
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图 9 CU 试验 的 孔 辽 水 压力 与 轴 癌 应 变 关 系 曲 线 
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Fig.9 Curves of the pore water pressure vs the axial strain in CU tests 


三 由 图 8 可 知 , 固 结 不 排水 试验 的 应 力 应 变 曲线 
相对 要 稍 小 一 些 。 这 是 因为 体积 的 不 变 会 影响 墙 体 
的 伺服 控制 ,从 而 影响 到 偏 应 力 的 大 小 。 由 图 6 和 
图 8 可 看 出 ,在 硬化 阶段 前 期 曲线 基本 都 呈 直线 状 
态 。 由 图 9 可 知 ,在 固 结 不 排水 模拟 试验 中 孔隙 水 
压 模 拟 曲 线 与 室内 试验 曲线 和 理论 计算 结果 也 较为 
一 致 ,孔隙 水 压 会 随 着 围 压 和 含 泥 量 的 变 大 而 变 大 ， 
并 会 出 现 负 的 孔隙 水 压 。 上 述 结果 表明 ,针对 所 研 
究 的 饱和 含 泥 沙土 建立 的 PFC3D 模型 是 正确 的 ,并 
可 作为 文献 [1] 理论 模型 验证 的 数值 试验 机 ,用 于 
理论 模型 的 真 三 轴 验 证 。 


3 ”饱和 砂 土 的 真 三 轴 CU 试验 结果 对 
比 及 中 间 主 应 力 影 响 系数 的 确定 


将 PFC3D 所 建立 的 模型 在 最 小 主 应 力 os 


200 kPa ,中 间 主 应 力 a. 2300 kPa 下 进行 了 固 结 不 排水 
试验 (CU 试验 ) 真 三 轴 数 值 模拟 ,将 所 得 PFC 模拟 结 
与 文献 [1 中 的 模型 理论 结果 进行 了 对 比 ,理论 模型 中 
间 主 应 力 影响 系数 b 值 分 别 取 0、.0. 25 ,0. 50,0. 75 和 
1.00。 图 10 和 图 11 分 别 为 真 三 轴 固 结 不 排水 试验 的 
轴 向 应 变 与 最 大 主 应 力 和 中 间 主 应 力 间 的 偏差 应 力 
ol -oa 以 及 轴 向 应 变 与 孔隙 水 压力 u 间 的 关系 曲线 。 
由 图 10 和 图 11 可 知 ,PFC3D 的 模拟 结果 也 能 
较 好 地 反映 饱和 含 泥 砂 十 的 真 三 轴 固 结 不 排水 加 载 
特征 ,其 变化 规律 与 文献 [1j 的 计算 结 末 有 较 好 的 
一 致 性 。 另 外 ,由 图 10 比较 可 知 ,PFC3D 的 模拟 结 
采 与 文献 [1j 的 理论 计算 结 采 在 中 间 主 应 力 影响 系 
数 b =0.50 时 ,符合 度 最 佳 ,其 误差 小 于 5% 。 由 图 
11 比较 可 知 , 当 含 泥 量 为 0% 时 ,PFC 模拟 曲线 与 
b 20.50 时 的 理论 计算 曲线 在 整体 上 均 符 合 较 好 ， 
其 误差 不 超过 5% ; 当 含 泥 量 为 5% 时 ,PFC3D 的 模 
拟 曲 线 与 5=0.50 时 的 理论 计算 曲线 在 轴 向 应 变 小 
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T 1296 时 都 符合 得 较 好 ,误差 小 于 3% , 轴 回 应 变 大 
于 12% 后 ,误差 有 所 扩大 ; 当 含 泥 量 为 1096 时 ， 
PFC3D 的 模拟 曲线 与 5=0.50 时 的 理论 计算 曲线 在 
加 载 前 期 ( 轴 疝 应 变 小 于 6% 左右 ) 符合 得 较 好 , 误 


差 小 于 6% , 轴 问 应 变 大 于 696 后 ,误差 有 所 扩大 。 
通过 综合 比较 图 10 和 图 11 可 知 , 对 于 南昌 地 铁 饱 
和 含 泥 砂 士 而 言 ,在 静 真 三 轴 条 件 下 ,文献 [1 理论 
模型 中 的 中 间 主 应 力 影响 系数 取 0. 50 较为 符合 。 
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Fig. 11 


4 三 循环 动 三 轴 和 动 真 三 轴 数 值 CU 试 
ME 


4.1 循环 动 三 轴 数 值 CU 试验 验证 


考虑 围 压 振动 频率 .振动 幅度 以 及 饱和 砂 土 含 
泥 量 等 因素 ,进行 了 PFC3D 循环 动 三 轴 CU 试验 模 
拟 , 并 横 回 对 比 了 文献 [11-2] 的 理论 计算 结 有 末 和 试验 
结果 。 在 围 压 为 200 kPa、400 kPa, 振动 频率 为 
0.5 Hz 1 Hz 和 2 Hz, 振动 幅度 为 30 kPa 40 kPa, 
50 kPa 和 70 kPa ,饱和 砂 土 含 泥 量 为 0% 5% 和 1096 
条 件 下 ,分 别 进 行 循环 动 三 轴 试 验 。 由 于 对 比 曲 线 
结果 类 似 , 限 于 篇 幅 仅 列 出 0% ,596 和 1096 的 饱和 
砂 士 含 泥 量 , 围 压 为 400 kPa ,振动 频率 为 1 Hz 振幅 
为 70 kPa 的 数值 模拟 结果 对 比 情况 。 取 试验 振动 次 
数 的 前 20 次 进行 分 析 。 
4.1.1 应 力 应 变 关系 

图 12 ~ 图 14 为 围 压 400 kPa .振动 频率 1 Hz 和 


2/%0 
(b) 含 泥 量 5% 


真 三 轴 固 结 不 排水 的 轴 辐 应变 与 孔 际 水 压力 间 的 关系 曲线 
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(c) 含 泥 量 10% 


Curves of the axial strain vs the pore water pressure in true triaxial CU tests 


振幅 70 kPa 时 , 含 泥 量 0% .5% 和 10% 的 砂 十 对 应 
的 应 力 应 变 关 系 。 

由 图 12 ~ 图 14 可 知 ,gq-el 关 系 均 表 现 为 随 着 循 
环 加 荷 次 数 的 增加 ,应 变 逐 步 紧 缩 . 动 前 切 模 量 逐 渐 
增 大 沛 回 圈 面 积 逐 渐 减 小 ,循环 稳定 时 割 线 模 量 、 
滞 回 圈 面 积 均 保持 稳定 。 究 其 原因 ,在 室内 试验 时 ， 
刚 开 始 进行 循环 人 荷载 时 砂 士 颗粒 分 布 、. 基 础 都 不 均 
匀 , 丛 载 传递 下 来 的 能 量 不 能 得 到 及 时 耗 散 ; 但 随 着 
谷 载 次 数 的 增加 ,颗粒 通过 移动 得 到 重新 排列 ,接触 
更 加 均匀 , 耗 能 系统 也 更 加 完善 ,于 是 便 出 现 了 前 期 
应 变 发 展 快 后 期 均匀 稳定 的 现象 。 相 同 条 件 下 , 饱 
和 人 砂 士 的 含 泥 量 越 高 ,相同 应 力 条 件 下 所 发 生 的 应 
变 会 更 大 一 些 。 
4.1.2 孔隙 水 压 的 发 展 

图 15 ~ 图 17 为 围 压 400 kPa 振动 频率 1 Hz 和 
振幅 70 kPa 时 , 含 泥 量 0% 5% 和 1096 KIE EXT 
的 孔 队 水 压 变 化 对 比 情况 。 
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(b) 理论 计算 结果 与 PFC3D 模 拟 结 果 对 比 
图 12 含 泥 量 为 0% 时 的 应 力 应 变 曲线 
Fig. 12 Curves of q-e, with 0% clayey content 
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(b) 理论 计算 结果 与 PFC3D 模 拟 结果 对 比 


图 13 含 泥 量 为 5% 时 的 应 力 应 变 曲线 


Fig.13 Curves of q-e, with 5% clayey content 
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图 14 会 泥 量 为 10% 的 应 力 应 变 曲线 
Fig.14 Curves of q-e, with 10% clayey content 
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(b) 理论 计算 结果 与 PFC3D 模 拟 结果 对 比 
图 15 含 泥 量 0% 时 孔 际 水 压 - 循 环 次 数 曲线 
Fig. 19 Curves of u-N with 0% clayey content 
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(b) 理论 计算 结果 与 PFC3D 模 拟 结果 对 比 


图 16 含 泥 量 5% 时 和 孔 际 水 压 -循环 次 数 曲线 
Fig. 16 Curves of u-N with 5% clayey content 
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(b) 理论 计算 结果 与 PFC3D 模 拟 结果 对 比 


图 17 会 泥 量 10% 时 孔 际 水 压 -循环 次 数 曲 线 


Fig. 17 Curves of u-N with 10% clayey content 


e hE 15 ~ 图 17 可 知 ,孔隙 水 压力 4 与 循环 次 数 

的 数值 模拟 结果 与 室内 试验 和 理论 计算 结果 在 振 
由 ES 天 小 和 发 展 趋势 上 都 比较 吻合 ,数值 模拟 的 孔 阶 
水 妖 相 对 值 比 室内 试验 大 , 比 理论 计算 小 。 孔 际 水 
压 均 在 首次 加 载 时 快速 升 高 ,随后 在 循环 荷载 的 作 
用 下 循环 往复 逐步 上 升 ,如 应 力 -应 变 关系 一 样 在 前 
期 发 展 较 快 , 当 达 到 循环 稳定 即 应 变 发 展 较 慢 时 , 孔 
隙 水 压 也 逐步 稳定 下 来 。 饱 和 砂 土 的 含 泥 量 越 大 ， 
所 模拟 出 的 孔隙 水 压力 也 会 更 大 一 些 。 这 一 定 程度 
上 是 因为 饱和 砂 土 中 孔隙 被 泥 颗 粒 填 堵 住 ,导致 孔 
WBUKIEZI ET 


4.2 动 真 三 轴 数 值 CU 试验 验证 


真 三 轴 模 型 采用 六 面体 模型 ,通过 用 Fish 语言 
代码 编程 实现 如 下 加 载 :在 相应 排水 条 件 下 ,对 
PFC3D 数值 模型 的 3 个 方向 同时 加 载 至 目标 最 小 主 
应 力 ;保持 一 个 方 癌 的 目标 最 小 主 应 力 大 小 不 变 , 在 
态 外 两 个 方 回 同时 施加 应 力 , 直 至 达到 目标 中 间 主 


应 力 为 止 ;等 待 固 结 完成 后 ,在 其 中 一 个 方向 进行 多 
次 加 载 ,保证 数值 试 样 不 会 发 生 突变 ,直至 数值 试 样 
前 切 破坏 。 其 余 设 置 和 常规 三 轴 模 拟 试验 类 似 。 

由 于 本 研究 第 3 节 已 经 表明 ， 静 真 三 轴 条 件 
下 ,文献 [1] 理论 模型 中 的 中 间 主 应 力 影响 系数 
取 0. 50 较为 符合 。 为 了 验证 动 真 三 轴 条 件 下 文献 
[1] 理论 模型 的 正确 性 ， 将 本 研究 所 建 PFC3D f 
型 作为 数值 试验 机 ， 将 其 模拟 结果 与 理论 模型 计算 
结果 进行 对 比 。 在 文献 [1] 理论 模型 中 同样 取 
b =0. 50 ,分 别 计算 含 泥 量 分 别 为 0% 、5% 和 1096 
饱和 砂 土 在 振动 幅 值 为 70 kPa、 振 动 频率 为 1 Hz, 
中 间 主 应 力 分 别 为 250 kPa 和 300 kPa， 最 小 主 应 力 
分 别 为 150 kPa 和 200 kPa 条 件 下 的 固 结 不 排水 理 
论 结果 ， 并 与 PFC3D 模拟 结果 作 比 较 。 限 于 篇 幅 
且 结 果 类 似 ， 本 研究 只 列 出 中 间 和 最 小 主 应 力 分 别 
为 300 kPa 和 200 kPa 时 不 同 含 泥 量 下 动 真 三 轴 应 
力 - 应 变 和 和 孔 辽 水 压力 对 比 结果 ， 如 图 18 ~ 图 20 
所 示 。 
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图 18 不 同 含 泥 量 下 的 应 力 -应 变 关 系 对 比 


Fig.18 Comparisons of the stress vs the strain under different clayey contents 
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图 20 理论 计算 值 与 PFC 模拟 值 对 比 


Fig.20 Comparisons between theoretical data and PFC simulation data 


由 图 18 可 知 :PFC3D 数值 模型 模拟 结果 与 常规 ”面积 减 小 并 趋 于 稳定 状态 ;不 同 含 泥 量 的 砂 士 在 循 
三 轴 试 验 结 末 类 似 , 均 表现 为 在 加 载 前 期 产生 较 大 PMISEBE HH. P ,其 应 变 的 发 展 具 有 很 大 的 区 别 , 含 泥 
RE IRR ED i ARR; EWR ” 量 越 多 ,应 变 发 展 速 度 越 快 .动静 切 模 量 更 小 、 沛 回 
数 的 增加 ,应 变 过 步 紧缩 、 动 亲切 模 量 增 大 、 沛 回 圈 疾 面 积 更 大 。 这 与 常规 循环 动 三 轴 的 表现 相同 , 较 
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小 刚度 的 泥 颗 粒 可 以 软化 整个 试 样 。 

由 图 19 可 知 ,PFC3D 的 动 真 三 轴 数 值 模拟 结 
与 理论 计算 结果 的 孔 际 水 压 发 展 趋势 和 变化 幅 值 上 
都 比较 吻合 , 均 在 加 载 初期 迅速 增长 , 当 塑 性 变形 稳 
定 后 孔 际 水 压 的 发 展 也 趋 于 稳定 。 相 同 条 件 下 数值 
模拟 结果 的 孔 际 水 压 的 绝对 值 比 理论 结果 偏 小 ,但 
基本 接近 ;不 同 含 泥 量 的 砂 土 在 相同 循环 加 载 条 件 
下 了 筷 际 水 压 的 变化 在 绝对 值 和 振动 幅 值 上 都 相差 较 
大 , 含 泥 量 越 多 孔隙 水 压 的 振动 幅 值 和 绝对 值 都 
更 大 。 

由 图 20 进一步 比较 可 知 : 含 泥 量 为 0% 时 ,理论 
计算 结果 比 模拟 结果 要 大 4.796 左右 ; 含 泥 量 为 5% 
时 ,理论 计算 结果 比 模拟 结果 要 大 6% 左右 ; 含 泥 量 
为 10% 时 ,理论 计算 结果 比 模拟 结果 要 大 996 左右 。 
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对 南昌 地 区 饱和 砂 土 进行 PFC3D 数值 试验 验 
证 ; 效 值 试验 为 单调 人 答 载 作用 下 不 同 含 泥 量 饱和 砂 
凸 外 庄 柔 性 膜 的 常规 三 轴 CD .CU 试验 和 真 三 轴 CU 
试验 ,循环 荷载 作用 下 动 三 轴 和 动 真 三 轴 CU 试验 。 
上 基体 结论 如 下 。 

1) 常规 三 轴 CD 和 CU 试验 的 三 方 对 比 结果 可 
知 5 在 固 结 排水 试验 中 ,在 相同 应 力 条 件 下 ,PFC3D 
数 箱 模拟 结果 相 比 其 他 两 个 结果 要 更 大 一 些 , 在 
纪 极 排水 条 件 下 则 相反 。PFC3D 常规 三 轴 数 值 模 
所 SK 罕 内 试验 和 理论 应 力 -应 变 、 孔 际 水 压力 关系 较 
E 。 

2) 真 三 轴 CU 试验 数值 模拟 曲线 与 理论 计算 曲 
线 的 对 比 结果 表明 ,整体 上 PFC3D 模拟 与 中 间 主 应 
力 影 响 系数 上 =0.50 时 的 理论 曲线 吻合 均 较 好 。 

3) 动 三 轴 CU 试验 结果 表明 ， 不 同 含 泥 量 饱 和 
砂 土 的 PFC3D 数值 模拟 、 室 内 试验 和 理论 计算 结 
果 的 应 力 应 变 曲线 均 较 为 一 致 ， 表 明了 本 研究 的 饱 
和 含 泥 砂 土 数值 模型 的 正确 性 。 循 环 动 三 轴 试 验 
中 ， 其 他 条 件 不 变 时 ， 饱 和 含 泥 砂 土 的 应 变 发 展 速 
度 会 随 着 含 泥 量 的 增加 而 不 断 变 快 ， 同 时 孔隙 水 压 
力也 越 大 。 

4) 动 真 三 轴 CU 试验 结果 表明 ,在 不 同 含 泥 量 
下 ,饱和 砂 士 PFC3D 数值 模拟 与 理论 结果 间 的 应 
力 -应 变 和 和 孔 队 水 压力 关系 比较 吻合 ,再 次 证 明 
PFC3D 模型 在 研究 砂 土 力学 特性 时 具有 较 好 的 适用 
性 。 饱 和 砂 士 的 含 泥 量 越 多 ,应变 发 展 速度 越 快 , 动 


S, UTE E S 7] , 沛 回 轿 太 才 、 扎 辽 水 压 的 振动 幅 值 和 
绝对 值 都 更 大 。 随 着 集 载 次 数 的 增加 ,应 变 逐 步 紧 
4i BIHER Te HL CK r Ind SLE T RN EG T E 
状态 。 
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